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Abstract 


A lattice theory is developed for composite materials in which 
the matrix and inclusions of the composite are idealized as a gridwork 
of thermoelastic bars attached to small rigid masses. This model 
generalizes \the spring and mass model prominent in the classical 
lattice theory of solids. Replacement of the springsby elastic 
bars provides extra internal degrees of freedom beyond the central 
force model of the lattice theory. The lattice representation leads 
to a system of finite difference equations of motion which can be used 
to solve dynamic and thermal stress problems numerically. Equations 
for stress and couple stress are obtained, which are identical in form 
to those of the linear theory of micropolar thermoelasticity. 



Introduction 


Consider a matrix containing inclusions which are of nearly equal 
size and which are spaced relatively uniformly (spacing L) , as shown 
in Figure la. Such a composite material can be represented by the 
lattice model shown in Figure lb, in which the matrix and inclusions 
have been replaced by a gridwork of massless, flexible bars attached 
to small, rigid bodies. The model so introduced represents a 
generalization of the mass-spring model well known in the theory of 
crystal lattice. By introducing elastic bars to represent the 
interparticle forces one provides extra internal degrees of freedom for 
non-central forces and couples. The model is therefore expected to 
encompass certain physical phenomena not present in the classical theory. 
Indeed this turns out to be the case. 

In order that the lattice model represents the actual composite 

material in some optimum sense, "effective" values must be selected for 

the masses and the mass moments of inertia of the bodies, the stiffness 

and the thermal expansion coefficients of the bars, etc. For example, 

the effective mass of each rigid body might be taken as the sum of the 

3 

mass of a single inclusion and the mass of a volume L of the matrix. 
Similarly, effective stiffness coefficients for each bar could be found 
by considering the volume of material between each inclusion and the 
elastic moduli of the composite (determined, for example, using the classic 
rule of mixtures). It is not, however, the purpose of this paper to 
discuss the procedure for selecting the parameters of the lattice model. 
Rather, the aim here is to develop the equations which govern the deformation 
of such materials. 
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EQUATIONS OF MOTION 


A repeating unit of a two-dimensional lattice model is shown in 
Figure 2. Interactions between nearest and next nearest neighbors are 
transmitted by means of thermoelastic bars. The bars are capable of 
supporting axial, shear and bending loads, as well as transmitting heat. 
The equations of motion for the rigid body at point (£,m) are 


».&,m & , m / i nv 

mu = F (a = 1,2) 


J = M 3 *’ m 


( 1 ) 


where u and are the translations and rotation of the lattice 

a Y 3 ^ ^ 

pbint (&,m), and F^ * and * refer to the resultant forces and moment 
exerted on the particle by the adjoining beam elements. The constants 
m and J denote the mass and mass moment of inertia of the body. These 
forces and moment result from both temperature changes 0 in the beams and 
from relative displacements (translations and rotations) of the end points 
of the beams. The forces which each element exerts on mass m may be 
calculated* from the stiffness-matrix and the thermal expansion properties 
of that particular element. It is assumed that all bars connecting nearest 
lattice points have identical properties, and likewise all bars connecting 
next-nearest points are identical. Summing the contributions of all eight 
bars gives the resultant generalized forces M^ ,m . Substitution 

of these expressions into (1) yields 


m u. 


£,m 


La^ (u 1 A+1, “-2u 1 A,m + 

+ Lb 1 (u 1 Jl,m+1 -2u 1 2 ' ,m + u^’” 1-1 ) 


*The details of computing the forces which a typical bar exerts on the mass 
are given in the Appendix. 
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I 


m u, 



. L , . , w &+l,m+l S,+l,m-l , £-1 ,m+l £-l,m-l . £,m. 

+ j (a 2 + b 2 )( Ul + u l + U 1 + U 1 - 4u l > 


, L , , N , £+1 ,m+l . £+l,m— 1 £— l,m+l £-l,m-lx 

+ 2 (a 2 V (u 2 " u 2 “ u 2 + u 2 > 


+ r b ! ^ 3 £,m+1 - ^ 3 " ,m “ 1 ) 


+ y- b 2 (*3 A+1 » mfl - ^ 3 *' +1 ’ m ' 1 + «j> 3 £ “ 1 ’ m+1 - <p 3 l ' 1,m 'b 


~ 2~ S 1 (9 Z+1 ’ m ~ 0 * 1,m ) 


_ g (Q^ +1 > m+1 + 0 A+l,m-l _ 0 £-l,m+l _ 0 A-l,m-lj 


£,m T / £,m+l 0 £,m £,m-lx 

? = La x (u 2 - 2 u 2 * + u 2 ) 


, / £+l,m 0 £ ,m . £-l,nu 

+ Lb 1 (u 2 * - 2 u 2 * + u 2 ) 


, L , , w £+l,m+l £+l,m-l £-l,m+l £-l,m-l . £,niv 

+ J ( a 2 + b 2 ) (u 2 + u 2 + u 2 + u 2 “ 4u 2 > 


. L , , w £+l,m+l £+l,m-l £-l,m+l £-1 ,m-lv 

+ 2 (a 2 " V (U 1 " U 1 - U 1 + U 1 > 


L 2 , / , £+l,m , £-l,m N 

— b x (<f> 3 “ *3 ) 


L 2 . ,, S,+l,m+l , , A+l,m-l , i-1 ,m+l , 

— b 2 (^ 3 + *3 - <t > 3 “ *3 > 


- £ 8l te 1 -”* 1 - e*’”- 1 ) 


,£+1 ,m+l £+1 , m- 1 £-l,m+l £-1 ,m-l> 

— g 2 (0 -0 +0 -0 ’ ) 


...( 2 ) 


" £,m L 2 , r . Jl+l,m £-l,nu , £,m+l 

J<f> 3 = — [ (u 2 - u 2 )-(«•!_ - u i )] 


+ L 2 r , i,+l,m+l £+l,m-l «,-l,m+l 

— h 2 [ (U 2 +U 2 - u 2 - u 2 ) 


/ £+l,m+l £+l,m-l , £-l,m+l £-l,m-lx 1 

- ( u i - u x + u - u x )] 


2 d l + + 4’3*’ mfl + <P 3 Z ’ m 1 + 8<}> 3 A ’ m ] 


L j r , £+1,111+1 . £+l,m-l , £-l,m+l , £— l,m-l Q , £,m-, 

- J d 2 *-^3 + ^3 + ♦3 + 93 + 8< p 3 J 


Here the stiffness coefficients , 


are defined as 


_ e l a l 

~ L-Lg 


12E I 

b = — £-£■ 
C - L*L_3 
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4E I 


<U E 


L-L 


L*L 


g 

H ~ L-L c 


£ = 1,2 (no suimnation) 


(3) 
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where E^, A^, I^, a^, and denote the Young’s modulus, cross-sectional 
area, moment of inertia, coefficient of thermal expansion, and the length 
of a bar t es P ect ively • A subscript £ = 1 designates bars connecting 
nearest lattice points; £ =* 2 refers to elements connecting next-nearest 
points. 

The equations of motion (2) can be written as differential equations 
by first expanding the displacements of neighboring lattice points in 
Taylor series about the point (Ji,m); i.e., 


£+l,m f , T » . « 2 \ 

u = (1 + L3-+ 3 n + • • . ) u 

a 1 21 1 a 


£,m 


u 


&-l,m 


= (1 - Wl +JT\ 2 - ...) u a 


l ,m 


(A) 




= (1 + L3^ + h^2 + 2T + L 2 3^3 2 + 2 T ^2 2 + ***^ U a 


ii » 2 




u 11+1,111 1 = (1 + L3_ - L3 + 3 2 - L 2 a.a, +^3, 2 + ...) u 

a 1 2 - 2! 1 1 2 2! 2 a 


£,m 


u = (1 - L 3 - L3, + ~ 3. 2 + L 2 3.3, + 77 3, 2 - ...) u 11,111 

a 1 2 2 ! 1 122! 2 a 

where ^ and ^ respectively represent the partial differential operators 
3/ 3x^ and 3/3x 2 . 

Substituting (4) along with similar expansions for the rotations 
into (2), and retaining only second derivatives, gives 


m 


. 3 U 1 " (a l + a 2 + b 2 ) u l,ll + (a 2 + b l + b 2 ) u l,22 


+ (2a 2 - 2t> 2 ) u 2 ^ 12 + (b 1 + 2t> 2 ) <f> 3 ^ 2 




(5) 
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m 

L 3 


(a l + a 2 + V U 2,22 


+ (2a 2 - 2b 2 ) u 1>12 


- (g x + 2g 2 )0, 2 


+ (a 2 +b 1 + b 2 ) u 2> 
’ (b l + 2b 2>*3, 


— (j> = (b + 2b )(u„ . - u „) 

L 3 3 1 2 2,1 1,2 


( 2 + d 2 ^ ^ 3,11 + ^ 3 , 22 ^ ( 2 b 1 + 4 b 2 ) i () 3 


.(5) 


Equations (2) and (5) provide a set of difference and differential 
equations, respectively, which govern the motion of the two-dimensional 
lattice model under consideration. A solution to these equations which 
satisfies prescribed boundary conditions yields the displacements u^ 
and the rotation <f>^. One can then compute the forces and moments in the 
lattice bars using the appropriate stiffness coefficients (given in the 
Appendix). Finally, stresses and couple stresses can be found by 
averaging the forces and moments over the faces of a small cube surrounding 
a lattice point. 


STRESSES AND COUPLE STRESSES 


To obtain the general constitutive equation for the material 
represented by this model we isolate the cube surrounding mass (£,m), 
as shown in Fig. 3a. This element is subject to the concentrated forces 
and moments transmitted by those bars which are intersected by the cube. 
These generalized forces are expressed in terms of the translations and 
rotations at the end-points of the beams, as described in the Appendix. 
Averaging the forces over the faces of the cube, and then making use of 
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(4) where only first derivatives are retained, gives the stress components 
(see Fig. 3 


fc ll = (a l + a 2 + V u l,l + (a 2 " b 2 ) u 2,2 • (8 1 + 2g 2 )0 

fc 22 = ^ a l + a 2 + b 2^ u 2,2 + ^ a 2 " b 2^ u l,l “ ^ 8 1 + 2g 2^ 0 

t 12 = (a 2 _ U l,2 + (a 2 + b l + b 2^ u 2 ,1 “ (b l + 2b 2' > *3 

fc 21 = (a 2 ‘ b 2 ) U 2,l + (a 2 + b l + V U l,2 + (b l + 2b 2 )<fl 3 


and the couple-stress components 
d l d 2 

m 13 ( 4 + 2 )lj> 3,l 


d l d 2 

m 23 ( 4 + 2 ^3,2 


It is of interest to compare equations (6) and (7) with the 
constitutive equations of Eringen's micropolar theory [1,2]. The stresses 
and couple stresses in a linear micropolar anisotropic thermoelastic 
solid are given by 


A £kmn £ mn + 

~~ C £kmn ^m,n 

where 

A £kmn ~~ A mn£k’ C £kmn C mn&k 


( 8 ) 


(9) 
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and 


“ e k& E kS,m ^m hP 

3 2 i(u k,£ + u A,k> 

1 

^m " 2 6 mnp U p,n 


( 10 ) 


In equations (8) and (10) <j>^ is the "microrotation", i.e., the average 
rotation of "microelements" about the center of mass of a deformed 
"macrovolume element". For the present model we can interpret the small 
rigid body at lattice point (£,m) as a microelement, and the surrounding 
cube as a macrovolume element. 

In view of (9) the number of coefficients (or C rt1 ) is 45. 

J&kmn iokmn 

For an orthotropic material (one having 3 orthogonal planes of symmetry) 
equations (8) become, in the two-dimensional case 


hi = hill hi + A 1122 e 22 + B ll 0 

fc 22 = A 1122 6 ll + A 2222 s 22 + B 22 0 

fc 12 “ A 2112®12. + A 2121 e 21 

hi = A 1212 6 12 + A 2112 e 21 

and 

m 13 = C 3131+3,l 

m 23 = C 3232+3,2 


( 11 ) 


( 12 ) 


The lattice model under consideration is symmetric also with respect to 
the two orthogonal 45° planes which bisect the planes of symmetry x^ - 
and x 2 - Xg 
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(see Fig. 2). Accounting for this symmetry and using (10), equations 
(11) reduce to 


fc ll = A 1111 U 1,1 + A 1122 U 2,2 + B ll 0 

fc 22 = A 1122 U 1,1 + A 1111 U 2,2 + B ll 0 

fc 12 = A 2112 u l ,2 + A 1212 U 2, 1 “ (A 1212 “ A 2112^ *3 

fc 21 = A 1212 u 1,2 + A 2112 u 2 , 1 + (A 1212 " k 2112)*3 
and equations (12) become 


m 13 C 3131^3,l 


(14) 


m 23 C 3131*3,l 

A comparison of (6), (7) with (13), (14) shows that the constitutive equations 
for the lattice model and the micropolar theory are identical providing 


A llll = a l + a 2 + b 2 

A 1122 = a 2 " b 2 

A 2112 = a 2 " b 2 

A 1212 = a 2 + b l + b 2 

B 11 = “ ( g i + 2g 2^ 
d l d 2 

n _ _i _i 

u 3131 42 


(15) 
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Hence, the material constants in the micropolar theory can be evaluated 
in terms of the stiffness and thermal expansion properties of the lattice 
elements. Considerable insight into the relative magnitudes of the 
material constants is thus provided. The approach also suggests a 
method for investigating how the elastic properties of a micropolar 
solid would be affected by changes in the geometrical and material 
properties of the constituents. In effect, it permits the construction 
of a material with a desired property, within the context of micropolar 
theory. 

The equations of motion for the element of mass (i.e., the cube) 
are given by the first and second laws of motion of Cauchy which express 
the local balance of momenta. 

Restricted to the two-dimensional case, in the absence of body forces 
and couples these are 


and 


t 3a,6 


pii 

a 


(ct ,3 = 1,2) 


(16) 


m 33,3 + e 3By t 6Y 


Pj < t>3 


(17) 


Substituting ( 6 ) into (16) gives equations (5^ providing 


P 


_ m 

= T? 


(18) 


Hence the displacement equations of motion for the lattice point 
(£,m) are identical to those for the cube element. 

Substituting ( 6 ) and (7) into (18) gives 


Pj 4> 3 = (b 1 + 2 b 2 )(u 2>1 - u 1>2 ) 

+ (^ + 2^)^ 3>1 i + ( f , 3 > 22 ) 


— (2b^ + 4t> 2 ) 3 


( 19 ) 
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Letting 


Pj =~iJ (20) 

we note that equation (19) then differs from (5^) only in the 
coefficient of the term involving (^ ^ 4- 22 } • This difference 

is due to the fact that the concentrated couples acting on the lattice 
point (£,m) are different from those acting on the cube element. Note 
that the cube is intersected by some lattice bars (e.g. bars (£,mfl - 
£4*1, m), (£4-1, m - £,m-l), etc.) which do not intersect the lattice 
point. These bars therefore transmit couples to the cube, but not to 
the lattice point. In addition there are couples acting on the cube 
which are associated with the transfer of shear forces from the lattice 
point to the faces of the cube. 
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APPENDIX : INTERACTION FORCES IN THE LATTICE MODEL 


The nature of the interaction forces transmitted by a thermoelastic 
beam in the lattice model is now examined. Each element is assumed to 
be a straight bar of uniform cross-sectional area A, flexural rigidity 
El, length L and coefficient of thermal expansion a. An element is 
capable of resisting an axial force, shear and flexure. The location and 
direction of the generalized forces (Q^) and the corresponding generalized 
displacements (q^) which act on a typical element are shown in Fig. 4. 

The forces are linearly related to the displacements and to the 
average temperature change A0 in the beam. The temperature is assumed to 
be constant over any cross-section, and to vary linearly along the axis 
of the beam. In this case the forces Q_^ are given by* 

{Q} = [S] {q} + [S t ]A0 (A. 1) 


where the stiffness matrix [S] is 



*A rectangular matrix is denoted by [ ] ; a column matrix by { }. 


(A. 2) 
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and the forces resulting from thermal expansion of the bars are 
determined using 


[s T ] 


aEA 

0 

0 

aEA 


0 

0 


(A. 3) 


We now consider a particular beam element, for example the one connecting 
lattice points (&,m) and (£ + 1, m + 1) as shown in Figure 5, P^, 

and Pg denote the generalized forces which act on mass (£,m), and 
P^, P^ and P^ denote those acting on particle (& + 1, m + 1) . These 
forces may be expressed in terms of the equal and opposite forces which 
the particle exerts on the beam 


{P} = [A] {Q> 


where from geometry 


1 

f~i 

l 

n 

0 


0 


0 


0 


1 

n 



0 


0 


0 


0 


0 


0 


- 1 


0 


0 


0 


(A. 4) 


0 

0 

0 

0 

0 

0 

0 

0 

0 

- _I_ 
~ fl 

1 

fl 

0 

1 

~ fi 

1 

" fl 

0 


0 0 - 1 
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The generalized displacements (assumed infinitesimal) may be 
expressed in terms of the displacements of the masses (£,m) and 
(£ + 1, m + 1) as 


/ i r-n i &,m a &+1 ,m+l £,+l ? m+l , &-KL,m-HL 

{q} = [BH^ * , u 2 * , <fr 3 * •, u x , u 2 , * 3 i 


where 


[B] = 


1 

n 

1 

n. 


1 

/~2 

1 

/~2 


0 

1 

n 

1 

n 

0 


1 

n. 

1 

n 


(A. 5) 


(A.5‘) 


The average temperature change in the beam is 


AQ = 




(A. 6) 


where 


[C] = [1/2 1/2] ' 


(A. 6’) 
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Substituting (A. 2) through (A. 6) into (A.l) gives 


r _-i r Trir , 2,,m &+l s nri-l , 

{P} = [K]{u 1 , u 2 , ^3 , u i , u £ , $ 3 } 


+ [K T ]{0 a,m 5 0 S ' +1 ’ m+1 > 


(A. 7) 


where 


[K] = [A] [S] [B] 
[Kj] - [A][S t ][C] 


(A. 7’) 


In expanded motion, the interaction forces transmitted to mass 
(£,m) become 

t, L . , , w S,+l,m+l &,m N 

P 1 = 2 (a 2 +b 2 )(u l ~ U 1 > 


+ 2 (a 2~ b 2 )(u 2 “ u 2 ) 


+ T h 2 <*3 




" 2“ s 2 (0 fi ' +1,nrf ' 1 + 0> £,m ) 

P 2 = 2 (a 2 + b 2 } ( U 2 S ' +1,mfl " U 2 A,m) 


(A. 8) 


, L / , w £+l,m+l 

+ I (a 2" b 2 )(u l ~ U l } 


L 2 , , , £+l,nri-l &,bk 

- ~ b 2 (<}, 3 + *3 } 


Li g2 (e ft+ 1 .«n + 0 ^ m ) 
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L 2 , r , £,+1 ,mfl &,m, , Jt+l,m+l £,,nu. n 

P 3 = r b 2 I(u 2 -“2 >- (u l -"i » 


- 2 d 2 <* 3 ’' +1 ’ n+1 + 2 * 3 t,m ) 


where 


V 2 

2 L*L„ 


= 


12E 2 I 2 

h 0 = 

2 L-L 3 


d„ = 


4E 2 X 2 


2 L*L„ 


(A. 9) 


a 2 E 2 A 2 
g 2 " L*L„ 


The subscripts 2 in (A. 8) and (A. 9) have been introduced in order 
to indicate that these quantities refer to a beam element which connects 
next-nearest lattice points. A subscript 1 is used to designate nearest 
lattice points. 

Expressions for the interaction forces in the other bars are found 
in a similar manner. 
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b) Corresponding Lattice Model* 
Figure 1 , 
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a) Cube Surrounding Lattice Point 


b) Stresses and Couple Stresses. 


Figure 3 , 



Figure 5 . Beam Element Connecting Lattice Points (£,m) and 
(JH-1 ,m+l) . 
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